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In a recent paper by Giuliani and Rothman |16j . the problem of 
finding a lower bound on the radius i? of a charged sphere with mass 
M and charge Q < M is addressed. Such a bound is referred to 
as the critical stability radius. Equivalently, it can be formulated as 
the problem of finding an upper bound on M for given radius and 
charge. This problem has resulted in a number of papers in recent 
years but neither a transparent nor a general inequality similar to the 
case without charge, i.e., M < 4i?/9, has been found. In this paper we 
derive the surprisingly transparent inequality 



The inequality is shown to hold for any solution which satisfies p + 
"^Pt P, where p > and pT are the radial- and tangential pressures 
respectively and p > is the energy density. In addition we show that 
the inequality is sharp, in particular we show that sharpness is attained 
by infinitely thin shell solutions. 

1 Introduction 

Black holes for which the charge or angular momentum parameter equals 
the mass are called extremal black holes. They are very central in black hole 
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thermodynamics due to their vanishing surface gravity and they represent 
the absolute zero state of black hole physics. It is quite generally believed 
that extremal black holes are disallowed by nature but a proof is missing. 
One possibility to obtain an extremal black hole is to produce one from the 
collapse of an already extremal object. Previous mainly numerical studies 
([E]) H) have concluded that when Q < M collapse always takes place 
at a critical radius Rc outside the outer horizon, and as Q approaches M, 
this value approaches the horizon. This is similar to the non-charged case 
where the Buchdahl inequality implies that collapse will take place when 
R < 9M/4, i.e., Rc = 9M/4, cf. [lO]. In the charged case the critical value 
is expected to be smaller due to the Coulomb repulsion, and this is also 
shown to be the case below and in particular as Q — > M the stability radius 
does approach the outer horizon. For more information on the relation of 
this topic to extremal black holes and black hole thermodynamics we refer 
to [7], [16], [13] and [11] and the references therein. 

The problem of finding a similar bound as the classical Buchdahl bound 
for charged objects have resulted in several papers; some of these are ana- 
lytical, cf. [E], [H], [13], [H], [E] and [20], whereas others are numerical 
or use a mix of numerical and analytical arguments, cf. [7], [12], and [15] to 
mention some of them. We refer the reader to the sources for the details of 
these studies but in none of them a transparent bound has been obtained 
(except in very special cases), on the contrary they have been quite involved 
and implicit. Moreover, most of these studies rely on the assumptions made 
by Buchdal, i.e., the energy density is assumed to be non- increasing and the 
pressure to be isotropic. 

In this work we will show that 



given that q < r (which is a physically natural assumption, cf. the discussion 
below), and that p + 2pT < p, where p > and pt are the radial- and 
tangential pressures respectively and /) > is the energy density. Here we 
have used lower case letters mg,q and r to stress that the inequality holds 
anywhere inside the object. We refer to the equations ([3]) and Q below for 
the exact definitions of these quantities. To the best of our knowledge this 
bound has not appeared in the literature before. 

In the non-charged case a general proof of the Buchdahl inequality 
2m /r < 8/9, in the case when p + 2pT < p, was first given in [Ij. A 
completely different proof was then given by Stalker and Karageorgis [18j 
where also several other situations were considered, e.g. the isotropic case 




(1) 
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where p = px- The advantage of the method in [18j (which is related to 
the method by Bondi [8] which however is non-rigorous) compared to the 
method in is that it is shorter and that it is more flexible in the sense 
that other assumptions than p + 2pT < p can be treated. On the other 
hand the result in [TB] is weaker than the result in [I] in the sense that the 
latter method implies that the steady state that saturates the inequality is 
unique, it is an infinitely thin shell. Indeed, in [1] it is shown that given 
any steady state, the value of 2m/r for this state is strictly less than the 
value 2m /r of a state for which the matter has been slightly re-distributed 
and this monotonic property continues until an infinitely thin shell has been 
reached for which 2m /r = 8/9. The method in [H] also shows sharpness but 
only in the sense that there are steady states with 2m /r arbitrary close to 
8/9, leaving open the possibility that different kinds of steady states might 
share this feature. Moreover, since the assumption p + 2pT < /c is satisfied 
by solutions of the Einstein- Vlasov system it is natural to ask if there exist 
regular static solutions to the coupled system which can have 2m jr arbitrary 
close to 8/9. This question is given an affirmative answer in [2j, where in 
particular it is shown that arbitrary thin shells which are regular solutions 
of the spherically symmetric Einstein- Vlasov system do exist. On the con- 
trary, the matter quantities and the corresponding spacetimes constructed 
in [18] for showing sharpness cannot be realized by regular solutions of the 
Einstein- Vlasov system. The construction in [18] gives that a solution which 
nearly saturates the inequality 2m jr < 8/9 satisfies |? + 2pT = p, and in ad- 
dition pt and p are discontinuous. Neither of these two properties can be 
realized by regular solutions of the (massive) Einstein- Vlasov system. 

In the present work where we study charged objects we will adapt the 
method in [18] to show the inequality ([1]) and its sharpness. This again sup- 
ports the claim above that this method is very flexible. We have not been 
able to carry out the strategy in [1| in this case. If we have succeeded it would 
have given a more complete characterization, cf. the discussion above. How- 
ever, we do show in Theorem 2 below that an infinitely thin shell solution 
(with properties specified in the theorem) saturates the inequality, although 
we cannot show that no other steady states can saturate it as well. We 
also mention that in [5] a numerical study of the coupled Einstein-Maxwell- 
Vlasov system is carried out which supports that there are arbitrarily thin 
shell solutions for this system which saturate the inequality ([1]). 

The outline of the paper is as follows. In the next section the Einstein 
equations will be given and some basic quantities will be introduced. In 
section 3 the main results are stated and section 4 is devoted to the proofs. 
In the final section we discuss our inequality in view of the bound derived 
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in [16j for a constant energy density profile. 



2 The Einstein equations 



We follow closely the set up in [16] but here we also allow the pressure to be 
anisotropic, i.e., the radial pressure p and the tangential pressure px need 
not be equal. We assume throughout the paper that p, the energy density p, 
and the charge density are non-negative. We study spherically symmetric 
mass and charge distributions and we write the metric in the form 

where r > 0, 6 G [0, 7r],(/? S [0, 27r]. It is well-known that the Reissner- 
Nordstrom solution for the charged spherically symmetric case gives 

e-2A(r) = i_^ + 9l = e^Kr) (2) 

Here R is the outer radius of the sphere and Q is the total charge. This 
solution is a vacuum solution. The purpose of this work is to investigate the 
behaviour of A and p when the matter and charge densities are non-zero for 
r < R. Before writing down the Einstein equations let us introduce some 
quantities following [16]. Let 



g(r) =47r / e^^+^^^^^T^Yd?/, (3) 

and 







r 
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mi{r) = 4:11 7] pdrj, (4) 
JO 

where q{r) is the charge within the sphere with area radius r and mi{r) is 
the mass within this sphere. The subscript i is used to distinguish mj from 
the gravitational mass mg which is defined below. Let us also introduce the 
quantity 

JO ^ 

The Einstein equations for A and p now read (cf. [7J and [16j ) 

- + — ^=87rp+^, (5) 

and 

2- = -87rp+ — 6) 
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where the subscript r denotes differentiation with respect to r. Equation ([5]) 
can be written as 



l-87Tr'p-^, (7) 



so that 



dr 

-2X^^_ 2m, (r) _ F(r) 
By requiring that dS]) matches the exterior solution ^ at r = R gives 



g-2A ^ ^ _ _ lAlz. (8) 



2M Q2 1 2 



or 



1 

M = - (8W + ^)dv + 
^ JO ^ 



2i?' 

which defines the total gravitational mass M. In view of this relation we 
now define the gravitational mass nig within a given area radius r by 

/X Fir) q^(r) 
rugir) = m,{r) + + (9) 

In terms of the gravitational mass we thus get 



-2A(r) ^ ^ _ 2mg(r) g^(r) 



e-^^^^^ = 1 - ^^^^^^ + (10) 

Let us also write down the Tolman-Oppenheimer-Volkov equation which 
follows from the Einstein equations, cf. [7], but note that in our case p is 
allowed to be different from px which modifies the equation accordingly 

Pr = ^,+-{PT-P)-{P + P)e'' + 47rrp - 4) • (H) 

3 Set up and main results 

The problem of finding an upper bound on the total gravitational mass that 
a sphere of area radius R with total charge Q can hold, or equivalently, 
to find the smallest radius Rc, referred to as the critical stability radius, 
for which a physically acceptable solution of the Einstein equations can be 
found, is formulated in [IB] as follows: 

A physically acceptable solution should satisfy 

P > 0, p > and /i > -oo, (12) 
< Q < M, R> R+. (13) 
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Here i?+ = M + \J — is the outer horizon of a Reissner-Nordstrom 
black hole. The quantities and q should satisfy 

mg{R) = M, q{R) = Q, (14) 
Q < rng, rUg + ^Jrn!^--^ < r. (15) 

We see immediately that these relations imply that q/r < nig/r < 1. We 
will in addition assume that the following condition holds 

p + 2pT<p. (16) 

The condition ()16p is likely to be satisfied for most realistic matter models, cf. 
[9], and in particular it holds for Vlasov matter, cf. 0] for more information 
on this matter model. 

Remark. In [T] and [3] the following generalization of this condition was 
imposed, namely that 

p + 2pT < rip, for some > 0. (17) 

However, in contrast to the non-charged case where a bound on M is given 
by a simple formula depending on $7 the simplicity is completely lost in the 
charged cased except when Vt = 1. Now, the case = 1 should be considered 
as the principal case, cf. [9j, and in the non-charged case it is when $7=1 
that the classical bound 2m /r < 8/9 is recovered. 

We are now almost ready to state our main result but first we define 
what we mean by a regular solution of the spherically symmetric Einstein 
equations. We say that ^ := {fi, X, p,p,pT, j^) is a regular solution if the 
matter quantities p,p,PT and are bounded everywhere and except 
possibly at finitely many points, p has compact support and the equations 
©, dS]), dS]) and (fTT|) are satisfied (where the matter quantities are C^) and 
the constraints ()12p and (llSp are satisfied. 

Theorem 1 Let ^ be a regular solution of the Einstein equations and as- 
sume that il6\) holds. Then 

Moreover, the inequality is sharp in the subclass of regular solutions for 
which pt > 0. 
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Let us immediately make a consistency check so that (jlSp ensures that 
the stabihty radius is strictly outside the outer horizon. Thus we wish to 
show that the inequality ()18p implies that e"^^^^-* = 1 — + % > 0, or 
equivalently that 



An elementary computation shows that this is true as long as 




which always holds. 

The proof of Theorem [1] relies on the method in |18j for the non-charged 
case. In the introduction we discussed the strength of this method but also 
its shortages; the question of uniqueness of the steady state that saturates 
the inequality (jlSp is left open, and the constructed steady states that nearly 
saturate the inequality cannot be a solutions of the coupled Einstein- Vlasov 
system (these issues were answered in [1] and [2] respectively). Furthermore, 
it is not completely obvious from the construction in [18] that these solutions 
approach an infinitely thin shell. This point also carries over in our proof of 
sharpness in Theorem [1] and we therefore find it natural to include a proof 
of the fact that an infinitely thin shell does saturate (llSp , although we have 
not been able to adapt the strategy in [1] to show that no other steady state 
can have this property. Furthermore, the numerical study in [5] supports 
that the maximizer for the spherically symmetric Einstein- Vlasov-Maxwell 
system is an infinitely thin shell. 

We therefore investigate a sequence of regular shell solutions which ap- 
proach an infinitely thin shell and adapt the method in [3]. More precisely, 
let ^'fe be a sequence of regular solutions such that pk, {pT)k,jk ^^'^ Pk have 
support in Denote by the total gravitational mass and by Qk 

the total charge of the corresponding solution in the sequence and assume 
that Q := limfc^ooQfc) and M = lim^^oo -^fc exist and that sup^g^/r < 1. 
Furthermore, assume that r'^pkdr — > 0, and r^{2{pT)k — Pk) dr ^ 
as k oo. 




In view of inequality (llSp this holds if 
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Theorem 2 Assume that {^k}'kLi ^ sequence of regular solutions with 
support in [R^ , R] with the properties specified above and assume that 

lim ^ = 1. (19) 

Then 



,^ ^/R [r g2 

M = - h\/— + — . 20 

3 V 9 3i? ^ ' 

Remark. That sequences exist with these properties, in particular the prop- 
erty (jl9p . has been proved for the (non-charged) Einstein- Vlasov system, cf. 
[2] (and [6j for a numerical study). The investigation carried out in [5j also 
supports that such sequences exist for the Einstein- Vlasov-Maxwell system. 

4 Proofs 

Proof of Theorem [1], As described above our method of proof is an 
adaption of the method in [18] to the charged case. Let a regular solution 
be given and let us define 

F{r) -2 
1~ 2r' 

and let 



mx{r) = mi{r) ^ — = - (21) 



2mA o 2 

X = , y = onr p, z = 



Note that the conditions ()12p and (llSp imply that 

X < 1, y > 0, and z < 1. (22) 
Indeed, the two latter bounds are immediate and the former follows since 



T5]) gives that > rni — (r — m^)^ so that 



^^2mA^2m|_g2^2m3_ - (r - mgf ^ ^ 



Lemma 1 The variables {x, y, z) give rise to a parametric curve in [0, 1) x 
[0, oo) X [0, 1) and satisfy the equations 

Sirr^p = 2x + x- z, (24) 
Svrr^p = y, (25) 



2^ x + y - z ^ . _. _^ [x + y - z 



\2 



8TrrpT = — -x + y-i-z + —- — , (26) 

2(1 — X) 4(1 — X) 

where the dots denote derivatives with respect to [i := 21ogr. 
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Proof of Lemma [I].- The proof is a straightforward computation using 
the Einstein equations ([5]) and ([6]) and the Tolman-Oppenheimer-Volkov 
equation pT]) . 



□ 

Now let 

(3(l-x) + l + y-z)^ 

■w{x,y,z) = . 

1 — X 

Differentiating with respect to /3 gives 
A — 3x + y — z 



w 



(l-x)2 



[(3x - 2 + y - z)x + 2(1 - x)y + 2(1 - x)i] . (27) 



Now, using the expressions of the matter terms given in Lemma [T] the con- 
dition p + 2pT < p can be written 

(3x-2 + y-z)x + 2(l-x)(y-i) < (28) 

where 

a = 3x^ - 2x + (y - zf + 2{y - z). 
From ([27]) and ^ it now follows that 

w = \{3x - 2 + y - z)x + 2(1 - x)y + 2(1 - x)z] 

(1 — x)^ '- 

<- '^~^''^\;' a{x,y,z). (29) 
2(1 — x)^ 



Since 0<x<l, y > and < z < 1, it follows that w is decreasing 
whenever a > 0, which implies that 

w < mayiw{x,y, z), (30) 

where 

= {(x, y, z) : < X < 1, y > 0, < z < 1 and a(x, y, z) < 0}. 

To solve this optimization problem we introduce s = y — z and note that 
max^; w{x, y, z) = max^;' w{x, s) where 

= {x, s) : < X < 1, s > -1 and a(x, s) < 0}. 
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It is straightforward to conclude that there are no stationary points in the 
interior of E\ so the maximum is attained at the boundary dE' of E' . The 
Lagrange multipher method leads to the following system of equations 

(1 - x)(6(l + s)+ 4(3x - 1)) + 2(1 + sf = 0, (31) 
x(3x-2) + s(s + 2) = 0. (32) 

From (j32p we have that = —2s — x(3x — 2) which substituted into (j3ip 
results in the equation 

(x + s)(l-2;) =0. (33) 

If X = — s we get from (]32p that 4s'^ + 4s = so that either s = = x or 
s = — 1 = — X. In the latter case we get w{\,—l) = 0, and the former case 
gives w{Q, 0) = 16. We thus conclude that w < throughout the curve. 
Since p > it follows from the inequality it; < 16 that 

This is easily seen to be equivalent to 

(^-^)^<1^. (35) 
Taking the square root of both sides and rearranging leads to 



( /?^_ ^^_> + «!)( /I + <0, (36) 

W 9 3 V 9 3r^^^ ^ 3 V 9 3r^ " ^ ' 

Since the second bracket is always non-negative and vanishes only if mg 
q = we have 

- /r r 



which is the first claim. 

To show sharpness we will construct a spacetime such that the corre- 
sponding curve from Lemma[T]intersects a small neighbourhood of (xg, 0, Zg), 
where Zg < 1 is a given ratio g^/r^ and Xg is the corresponding value of x 
when equality holds in (f371) . i.e., 



_ 4 Zq 4/1 Zq 

9"y + 3V9^y 

We will construct such a spacetime by showing that there exists a curve 

X = x(r), y = y{T), z = z(t); r G [0, oo), 
which passes near (xg, 0, Zq) and in addition has the following properties 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

X 



Figure 1: The curves 70, 71, 7e and /i^ 



• (Al) is negative and locally integrable, 
. (A2) x(0)=y(0)=z(0)=0, 

• (A3) < x{t) < Xq, z{t) < Zq, 

• (A4) ^(t) = for all large enough r, x{t) — > and z{t) ^ as 
r — > cxD, 

• (A5) the curve is except for finitely many points. 

Below we will denote s = y — z as above and the curves {x{t), y{T), z{t)) G 
[0, 1) X [0, 00) X [0, 1), and (x(r), y(r), s(t)) E [0, 1) x [0, 00) x (-1, 00) wih be 
used interchangeably. Let us first see that if we have a curve which satisfies 
(Al)-(A5) a spacetime can be constructed. Indeed, let 

k(t) - dw 2(1 - xf 

a{x, s) dr 4 — 3x + s ' 

and observe that k is positive and locally integrable by (Al) and (A2). Next 
define 

P= Kdr, (39) 
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and 

r = (40) 
and define the metric coefficients by 

A = -^log(l-i) (41) 

It is straightforward to check that A and ^ solve the Einstein equations ([5]) 
and ([6]). The definition of k now imphes 



1 dw 4 — 3x + s 
'k'Jt ~ ~ 2(1 - x) 



where we recall that dots denote differentiation with respect to /? = 2 log r. 
Using ([27]) we thus have 

(3x + s-2)i; + 2(l = -^^y^, (44) 

which is equivalent to the relation p + 2pT = p in view of (j28p . We will now 
show that such a curve exists. Let us fix some small e > and define 

,):^»^-^-"^ + (45) 

1 — X 

Consider now the curve 7^ in the (x, s)-plane defined by 

w^{x, s) = (eVl + 3x + 4(1 - e)f. (46) 

Define the corresponding curve in by 

(x, y, z) = (x, max (0, s), max (0, — s)), (47) 

so that s = y — z. In figure 1 the curve 7^ is depicted together with the curves 
7o and 71. Note that 70 is the curve w{x, s) = 16 which passes through (0, 0) 
and (1, —1). The curve 71 is the curve a(x, s) = which also passes through 
(0,0) and (1,-1). The dotted line shows the line s = Sg := —Zq (for the 
choice Zq = 0.6) and the intersection of the curve 70 with this line is the 
point {xq,Sq). It is clear that for a sufficiently small e > the curve 7^ 
intersects an arbitrarily small neighbourhood of (xg,Sg). Let us denote the 
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point of intersection of 7^ and the line s = Sq hy (x^, Sg). Let us now define 
r := 7e + /le, where /i^ is the curve given by the equation 

— = — —, such that s{xl) = s„. (48) 
dx s + X ^ 

It is clear from the defining equation that £ {{x, s) : x > 0, s < 0, s + x > 
0} and that the solutions approach the point (0, 0) for all admissible starting 
points {Xq, Sq) (note that Xq + Sq > 0). The curve is depicted in figure 1. 

It remains to show that (Al)-(A5) are satisfied for the curve T and that 
p,p and pt are non- negative along the curve. It is obvious that T satisfies 
(A2)-(A5). To see that it satisfies (Al) we first consider the first part of the 
curve 7e and note that a > along 7^. This follows since 7^ lies above 71 
and Q = along 71 and 

da 

— = 2s + 2 > 0, for s > -1. 

OS 

Hence it is sufficient to show that dw/dr < to establish that 

1 dw 

-— <0. 

a dr 



We differentiate (jl6|) and obtain 

dwe 3e 3(1 — e)(l — x) + 1 + s dx 

dr ~ ^1 + 3x V'^-x dr' 

If we now differentiate (j45|) directly we get 



(49) 



dwe (3 — 3e)(l — x) + 1 + s r , , , .dx , .ds-, 

-Jf = TTTij^ - 3.)(1 - .) + 1 + + 2(1 -.)-]. 

(50) 

Comparing ([^9]) and ([50]) gives 

, /3e(l — x)^'^ , , dx , , 

Thus differentiating w along 7^, substituting for ds/dr using ([5T]). leads to 



dw _ 3e(4 — 3x + s) y/1 — x — \/l + 3x dx , , 

dr ~ 1 - X VI + 3x dr' 

Since dx/dr > along 7^ and since < x < 1 and s > — 1 we get that 

dw 
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It remains to show that a ^dw/dr is negative also along the curve h^. Here 
we have 

1 dw 3(1 — x) + 1 + s r ds-,dx 

" {x{3x - 2) + s{s + 2))(1 - xf ~ ^ ^ ^ ""^d^^d^ 
_ [3(1 -x) + l + s] dx 

{x + s){l-xy dr' ^ ' 

where we used (|48p for ds/dx. Since dx/dr < along the claim follows 
since x + s > along /ig. Hence, 

< 0, (54) 

a T 

along r. The local integrability of this expression follows by inspection of 
the formulas above since < x < x^ < 1. Thus condition (Al) holds along 
r. Finally we show that p,p and pT are non-negative along T. Since y > 
along r, cf. (j47p . it immediately follows that p >0. Since (j44p implies that 
p = p + 2pT we only need to show that pt ^ along T. First we consider 
the first part of F, i.e., the curve 7^. From Lemma[T]and ()5ip we have 

o 2 X + S . . ix + sY 

onr pt = —, rX + s — z -\ ; 

^ 2(1 -x) 4(1 -x) 

/3e^/l - X , 3 -.Idx . (x + s)^, , ^, 

= ( — : + 1 e) — :r + \-z + —, -r]- (55) 

^2VTT3^ 2 UdT ^ 4(1 -x)^ ' 

Since e is small the first term is positive since dx/dr > along 7^. The term 
in square brackets can also be seen to be positive. Indeed, along the part of 
7e where s > 0, z = and the claim is trivial so we focus on the part where 
s < 0. Here z = —s and we thus want to show that 

^^{x + sf_ 
4(1 - x) 



along the part of 7^ where s < 0. We evaluate the left hand side of (|56|) along 
7o and show that it is positive there and then we conclude by continuity that 
this statement also holds along 7^ for e small. Along 70 we have the relation 



4 + 3s 4/1 s 



A straightforward calculation now gives that along 70 



(x + sV ,A n 16 
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Hence pT > also along 7^ for a sufficiently small e. Along /i^ it holds by 
construction, cf. (j48p and Lemma ([1]), that /> = and thus p = pT = 0. This 
completes the proof of Theorem [TJ 



□ 



Proof of Theorem [2l Let us begin with a few general facts. Consider 
a regular solution where the matter quantities are supported in [0,i?]. We 
recall from section 2 the following consequence of the matching condition 



.-M.,M.)= A _H^£ + (58) 

so that e^'^^ = 1 for r > R. Let us now derive an explicit expression for 
The Einstein equation ([6]) can be written as 

/^^=(^ + 4vrr-p-|^ + ^)e2\ (59) 

so that 

y+4-p-|:3+^)e2^dr, (60) 

since /i ^ as r ^ cxd in view of ()58p . We will also need an explicit formula 
for Xr, and from ([5]) we have 

A. = (4.rpir) - ^ + |^ - |,)e-. (61) 
From the expressions of fir and we also obtain 

/•oo 

fi{r)+X{r) = - 4:TTr]{p + p)e^^dr], (62) 

Jr 

SO that in particular 

^ + A<0. (63) 

Now we derive our fundamental integral equation which is a consequence of 
the Tolman-Oppenheimer-Volkov equation. Let 

iP = (mg + inr^p - —)e^^^. 

Taking the derivative of V' with respect to r, a straightforward calculation 
using the Tolman-Oppenheimer-Volkov equation (llip results in the following 
equation 

(m„ + Attv^P - —)e^+^ = e^'^^UiTrj'^ip + p + 2pt) + %)dr]. (64) 
r Jq 1]^ 
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This equation must be satisfied by any spherically symmetric static solution 
of the Einstein-Maxwell system. 

Let us now consider our sequence of solutions. Since Pk{R) = 0, {mg)k{R) 
Mk, qk{R) = Qk and e^'=+^'=(i?) = 1 we get in view of ([M]) for r = R, 

^k-^= r e'^'^^^'^iiWipk+Pk + 2{pT)k) + 4)^^- (65) 

R Jr„ t 

Here we also used the fact that the matter is supported in [Rk, R\- We split 
the right hand side as follows 



[ e^'>'+^''{4TTri''{pk+Pk + 2{pT)k) + %) dTj 
= [ e^'^+^^{8Wpk + %dv 

rR 

+ / e^'>^+^'{47rr]\pk + 2{pT)k-Pk)d7]=:Sk + Tk. (66) 
JRk. 



'Rk 

By the mean value theorem we get that there is a ^ E [Rk, R] such that 



5fc = 2e'^'=«)e / e^'=(47rr?pfc + ^)dr/ (67) 
JRk 



=:Sl + Sl (68) 

Here we used equation ([^T|) for A,.. Now, since sup;, qk/r is strictly less than 
one we obtain a uniform bound on from the inequality ([T|), cf. the con- 
sistency check after the formulation of Theorem [TJ The same computation 
guarantees that {mi)k{r)/r + Fk{r)/2r < 1/2, thus it follows that 

R 

0<Sl< Clog (— ) ^ as A; ^ oo. 
Rk 

Since ^fc + < by ([63|) and since Pk > Pk + 2{pT)k > 2(j)T)k, it follows 
from the assumptions on the sequence that also 

Tfc — > as A; ^ oo. 

For the term Sl we get by (|58l) 



jRk di] ]j R 



R-^ > 
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Note here that Xk{Rk) = due to the support condition of the matter terms. 
In view of (j60p and the general bounds on nii/r and q/r it follows that 



g/^au ^ ^/ 1 _ _ + ^ as /c ^ oo, 



so that 



fc^oo " V i? i?2V Y R E? 

In conclusion, from ()65p we get in the limit k ^ oo, 



lini5^ = 2Wl-— + ^(1- Jl- — + ^). 



After some algebra this relation can be written as 



02 02 / Q2 

Squaring both sides one finds after some rearrangements 

(9M - + -^)(— - - 4Mi?(— - -^), 

so that 

(3M - = AMR. (70) 
it 

We have thus arrived at the same expression (with equality instead of in- 
equality) as (j35]) and we accordingly obtain 



VM = h a/ — + — , 

3 V 9 3ii' 



which completes the proof of Theorem [2j 

□ 



5 Final remarks 

In [18] several different conditions on the relation between p, p and px are 
investigated, e.g. the isotropic case where p = pT- We have not tried to 
consider other cases than p+px < /o in this work although we believe that it 
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can be done. We believe however that an equally transparent inequality as 
([1]) is unlikely to be found under other conditions than p + Pt ^ P, cf. |18j . 
However, the following comparison with the non-charged case is interesting. 
The original Buchdahl inequality [IQl was derived under the assumptions 
that p is non-increasing outwards and the pressure is isotropic and the steady 
state that saturates the inequality 2M/R < 8/9 within this class of solutions 
is the one with constant energy density for which the pressure is infinite at 
the center. It is quite remarkable that exactly the same inequality holds 
much more generally [1], as long as p + pT < p, and in particular that the 
steady state that saturates the inequality in this class is an infinitely thin 
shell which is drastically different from the constant energy density solution. 
One can now ask if there is a similar analogue in the charged case. 

In the work [16] by Giuliani and Rothman they find an explicit solution 
with constant energy density and constant charge density and they obtain 
for this solution an algebraic equation from which the values of the stability 
radius can be evaluated. It is in view of the discussion above therefore 
interesting to see whether these values are less, equal or greater than the 
values given by ([I]). In [16] the ratios R/M are displayed for different ratios 
oi Q/R (or more precisely for different ratios Q/M, but the corresponding 
ratios Q/M can be deduced). It turns out that the critical stability radius 
given by the relation 



are smaller than the corresponding ones found in [16], or alternatively, our 
relation admits a larger ratio M/R for a given ratio Q/M. 



I would like to thank the authors of [16] for their clearly written paper which 
got me interested in this topic. 
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